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Abstract 

We calculate the decay rate of light scalar mesons, using a diquark chiral effective theory, recently 
proposed to describe exotic hadrons. In the effective theory the light scalar mesons are postulated 
to be bound states of diquark and anti-diquark. We find our results are in good agreement with 

V^-"' experimental data. The axial couplings of diquarks with pions and kaons are found to be quite 

small and the perturbation is reliable. It shows that the diquark picture captures the correct 

^ \ physics of light scalar mesons. 
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I. INTRODUCTION 



Recent discovery of exotic baryons [l| revived the interest in hadron spectroscopy, es- 
pecially that of exotic hadrons. Among several models for exotic baryons, promising one 
is the diquark model 0, HI, which assumes a strong correlation of two quarks in the color 
antitriplet channel. In the Jaffe-Wilczek (JW) diquark model 0J, the exotic baryons are 
bound states of two scalar diquarks and one antiquark. The salient features of JW model 
are the occurrence of low-dimensional multiplets and small mass difference among the mem- 
bers in the multiplet, which were soon supported by the NA49 experiment 0]. Furthermore, 
the extremely narrow decay width of exotic baryons was naturally explained in the diquark 
picture of exotic baryons, where the decay process is suppressed by tunnelling among two 
diquarks [5||. 

As the existence of hadrons containing quarks more than three is logically possible, its 
absence in the low-lying hadronic spectroscopy has been a puzzle since the advent of the 
quark model p). In late seventies, Jaffe suggested that the positive-parity scalar mesons of 
mass less than 1 GeV might be bound states of multi quarks, forming a SU(3) nonet Q. 
In his scheme, the scalar nonet is postulated to be the non-exotic components of 4-quark 
states (Q 2 Q 2 ). However, the quark content of light scalar mesons has been controversial and 
still is |8[. Recently, it was argued that ex, the lightest member in the scalar nonet, is not a 
pure QQ state but mostly a four-quark state, using the unitarity argument in the large Nc 
analysis 0. (See also a recent paper for the diquark picture of scalar mesons (Io|.) 

In this paper, we assume that the scalar nonet is a bound sate of a diquark and an anti- 
diquark and calculate the decay width of scalar nonet to compare with the experimental 
data 11, using the diquark chiral effective theory (if. 



II. DIQUARK CHIRAL EFFECTIVE THEORY 

The success of quark model in hadron spectroscopy is understood in the framework of the 
chiral quark effective theory (xQET), derived by Georgi and Manohar [13]. The xQET is 
an effective theory of quantum chromodynamics (QCD) below the chiral symmetry breaking 
scale A xSB but above the confinement scale, Aq CD . The relevant degrees of freedom of the 
effective theory therefore contains quarks, gluons and the Nambu-Goldstone (NG) bosons 
of spontaneously broken SU(3)l x SU(3)r chiral symmetry. The chiral quark effective 
Lagrangian is given as 

C xQ = $(i$)+y)iP + g A $ 4 75^ - + l -f tr - ~ tr F^ v + • • • , (1) 

where = <9 M + ig s G^ is the QCD covariant derivative, 

V, = ~ + , Ap = \i (?d„£ - C^) (2) 

and / = F n = 93 MeV is the pion decay constant, 

S = e 2in/s = e 2 . (3) 

The Nambu-Goldstone fields n = ir a T a with tr (T a T b ) = l/2 5 afe . In xQET the higher-order 
terms or explicit chiral symmetry breaking terms are suppressed by powers of A x sb — 47r/. 
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The typical size of symmetry breaking terms is 



2 

T5-«0-l- ( 4 ) 

A XSB 

Therefore, they make only small contributions to the naive quark model. Furthermore, by 
matching with QCD, one finds the axial coupling and the strong coupling to be small as 
well, gA — 0.75 and a s ~ 0.28 so that the perturbation is reliable. This is why the naive 
quark model is so successful in describing non-exotic hadrons 

Similarly, the success of diquark picture for exotic hadrons may be explained by the 
diquark chiral effective theory (D^ET), proposed recently by Sohn, Zahed, and one of the 
authors . Since the diquarks have mass around 450 Me V 0, 13 , they fit naturally in the 
range where the xQET is applicable. Hence, in addition to the relevant degrees of freedom 
of xQET, the diquark chiral effective theory involves one more degree of freedom, namely a 
diquark field, defined as 

= J™ W J ^ k e^(xhMy), (5) 

where k is a mass scale for the diquark field, 7 m is the anomalous dimension of the diquark 
correlator and the Greek indices denote colors, while the Latin indices k = 1, 2, 3 denote 
flavors. i/; c = Cip T is a charge conjugated field of a quark field, if), an d C = 27270- (Here 
we consider only scalar diquarks, since tensor diquarks [3| do not contribute to light scalar 
mesons.) 

The diquark chiral effective Lagrangian is described by 

C = i|(23 M + *%)^| a -i^|^| 2 -(p^^*^^ 7b ^ + h.c.) 

+ h A Tr(cp* ( pAl)+ti A ip*ipTrAl + C mt (<p, ip, $ + C xQ , (6) 

where higher-order interactions of diquarks and quarks is denoted as £i nt , which may contain 
the (chirally) covariant derivatives. In this paper, we will take M = 420 MeV from the 
random instanton model 0| and g 2 = 3.03, the best value fitting the scalar nonet mass 0. 
The Yukawa coupling splits diquark masses, M u d = 450 MeV and M us = Md s = 630 MeV. 
And the binding energy between a diquark and an antidiquark in scalar nonet is B = 
280 MeV. We list these parameters in Table |U 



M g8 (MeV) M Md (MeV) B (MeV) g 2 
630 450 280 3.03 



TABLE I: The best-fit parameters in [jj. They are obtained from the random instanton model and 
the scalar nonet masses. 



III. DIQUARK- ANTIDIQUARK TO GOLDSTONE BOSONS 

In general, it is very hard to deal with bound states of strongly interacting particles, 
especially with broad resonances. However, we will attempt to calculate the decay width 
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FIG. 1: The annihilation process, <p + tp* — > tt + ir. The double line denotes diquarks and the 
dashed line denotes pions and kaons. 

of strongly interacting scalar nonet, utilizing the di quark chiral effective theory. We will 
follow the scheme of nonrelativistic positronium decay, assuming the diquark and NG boson 
couplings are not large. We will see later indeed the couplings are small and the naive 
diquark model of Jaffe and Wilczek is reliable. In the effective theory, the scalar mesons 
are bound states of diquark and anti-diquark. The decay process occurs as annihilation of 
diquark with anti-diquark into NG bosons (see Fig. [TJ) . The decay width is then given as 

r s _ Mnr = lim v 2 a(<p + V * ^7r + 7r)|^(0)| 2 (7) 

where v 2 is the velocity of anti-diquark in the rest frame of the target diquark and ^(0) is 
the IS wave function of the diquark-antidiquark bound state at the origin. The differential 
cross section for the annihilation process is then 

da = _ d$( k h) (g) 

^( Pl -p 2 )i-M?Mi 

where M. is the annihilation amplitude for <p* + <p — > 7r + n. The phase space is given as 

d 3 k d?k 

a$(p 1+P2 ; h, k 2 ) = 6\p 1+P2 -h- fa) (27r) B^ (27r) 3 2 a g4 - ( 9 ) 



Integrating over the phase space, we get 



S — >7T7V 



16tt M x M 2 1^2 ~ m 2 A + Jfi - m | 



|^(0)| 2 , (10) 



where Mi is anti-diquark mass, M 2 diquark mass, rriA the mass of pion A, tub the mass of 
pion B, and k\ (or — hi) is the final momentum of pions in the CM frame. 
To the order we are interested in 

M = Mtrcc + Mn oop . (11) 
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(a) (b) 
FIG. 2: The diquark-pion vertices: (a) vector coupling and (b) axial coupling. 

A. Tree contributions 

In D^ET, the leading order contribution to the amplitude in 1/A x sb expansion, shown 
in Fig|3 is given as 



7r A (k 1 ),n B (k 2 ) d ItpMMfa)) 



where 



A = Tr V" {<f*d^ - ipdtf*) + h A Tr + h A <P*<P Tr ^ 



(12) 



(13) 



Note that because of the parity invariance the leading coupling of (scalar) diquark with 
pions occurs at 1/A^ SB . For the manifest 577(3) covariance, we rewrite the incoming state 
in the basis of S77(3) octet and singlet: 



\ Va (jp 1 )^l(pz)) = 2Tg S c ( Pl ,p 2 )) + -5 ab S°( Pl ,p 2 )). 
Then, the tree amplitudes for octet and singlet become, respectively, 



M 



C-^AB 
tree 



4/ ; 



;{ki - k 2 ) ■ (pi -p 2 ) f 



ABC 



h, 



2f 



k\ ■ k 2 d 



ABC 



P 



(14) 

(15) 
(16) 



B. One-loop contributions 

The vector-current of pions (Fig. Et) contributes in the modified minimal subtraction 
scheme as 

-Mlloop = -^pf ABCT ^a P2 - h) ■ ( Pl - P2 ) h + ■ ■ ■] , (17) 

where the ellipsis denotes the higher-order terms and 

(xpi — yp 2 ) 2 — xp\ — yp\ + m 2 



h 



1 rl—x 

dx (4 — 3x + 3y) log 
o Jo 



(18) 
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in the SU(3) limit. The axial-current contribution (Fig. Eb) becomes 



Mi 



loop 



{h ■ k 2 ) \e alm e bon TiT* I 2 - {5 AB 5 ab - 2{T A T B ) ba }h + • • • ] (19) 



ArpB\ 



16n 2 f 2 



where 



rl rl—x rl—X—y 

1 2 = 6 / dx / dy dz log 

Jo Jo Jo 

rl rl—x rl—x — y 

1 3 = 2 dx dy dz log 

Jo Jo Jo 



(xpi + ypi — yk\ + zk 2 ) 2 — xpf — y(p\ — ki) 2 — zk\ + m 2 



(xki + yk\ + yk 2 + zpi) — xk\ — y{k\ + k 2 ) — zp\ + m 
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(a) 



(b) 



FIG. 3: One loop amplitudes: (a) vector coupling and (b) axial coupling. Solid lines denote quark 
fields. 



IV. TETRAQUARK SCALAR DECAY WIDTHS 

A. Decays into two nonstrange Nambu-Goldstone bosons 

In our constituent diquark picture, the light scalar mesons have a binding energy of 
280 MeV, which is not negligible in calculating the decay width. Because of the large binding 
energy, Eq. (fTU|) is far from the proper energy-momentum conservation. (The velocity v 2 
does not also go to zero. However, it factors out in the leading order.) To take into account 
of the binding energy in the energy-momentum conservation we modify Eq. (fTU|) as 



1 l-MI 2 I ifei I 



l*(0)| 2 , (20) 



\^MtM^P i+m 2 A + ^P i+m 2 B 

where M* = M X - BM l /(M 1 + M 2 ), M* = M 2 - BM 2 /(M 1 + M 2 ), and 

p _ {(M* + M*) 2 + m| - m 2 A } 2 



By a naive dimensional analysis we find 

|^(0)| 2 = c(MB) 3/2 (22) 

where c is a constant and M is the reduced diquark mass. From the analogy with the IS 
wavefunction of positronium, we take c—1/tt here. 
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The lightest scalar a is believed to have no or just little strangeness in it. So we assume 
the ideal mixing between a and fo- In our calculation, the (anti)diquark in o is nonstrange 
and the (anti)diquark in f includes one strange (anti)quark. We take = 0.75, the 
experimental value in the quark model k| 3], and set the quark mass m = 400 MeV [l6| . 



We then determine the couplings h A and h' A by fitting the experimental values for the decay 
widths of scalar mesons in the dominant modes. The results are shown in Table ITD [17]. We 
find our couplings are of order 1 for the decay widths in the range allowed by experiments. 
Since the expansion parameter is h\/Air or h'^/Air rather than the coupling itself, the higher- 
order corrections are quite suppressed and the perturbation is reliable. 



set 


h A 


h' A 


r CT ^(MeV) r /o _ 


.2. (MeV) 


r a „(MeV) 


1 


1.46 


-1.66 


317 


67 


83 


2 


1.20 


-1.38 


234 


71 


57 



TABLE II: Two fitting parameters and the decay widths of the light scalar mesons in dominant 
modes. 



Once we determine the axial couplings Ha and h' A , we calculate in the next section the decay 
width of fo in its subdominant mode to see if D%ET works. 



B. fo decay into two kaons 

Since the mass of two kaons is larger than fo mass, we cannot obtain the decay width with 
Eq. (f2T)|) . The broadening of fo makes the decay into two kaons possible. So the broadening 
should be considered in order to calculate the decay mode into two kaons. Keeping the basic 
picture, we introduce the Breit-Wigner form to revise Eq. (J201); 

r A B=N f Smax ds N rnsTs |^(0)| 2 (23) 

S ^ A * B J*™ 4tt 2 s + m 2 A + ^2 + m | (s - m|) 2 + m|T| 

with scalar mass m s and total decay width T S - y/s max is ms + Ts and y/s m i n is chosen 
to be the larger one between the threshold value for the decay mode and m$ — Ts- So 
y/s m i n = rriK + rax for f — > K + K mode. The wavefunction at the origin and the final 
momentum are modified as 



and 



|# s (0)| 2 w {(Mi + M 2 - yi)M} 3/2 /vr (24) 



k) = (s + ^b -O 2 _ m% (25) 
J 4s 



in (j23J- The normalization constant iV is adjusted by 



J_ r s ™«* 1 m s T s , , 

N J Smin aS 'tt (s - m|) 2 + m^l' 1 ' 
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Since /o — > 2ir is dominant among f decay modes, we will use Tf -+2n in Table ITT1 for 
in Eq. (j53J). Then we obtain 75 MeV with set 1 and 59 MeV with set 2 for Y^^k- 
Experimentally r f „-,•>. k / (T f + Tj ^2k) is 0.14 ~ 0.32 and the full decay width of fo is 
T = 40 ~ 100 MeV |l2j. We find our result is a few times larger than the experimental 
value, Tf ^2K = 5.6 ~ 32 MeV. Considering the experimental uncertainty and the crude 
approximation we made for the bound state wavefunction, our result is not too bad. 

V. CONCLUSION 

We calculate the decay width of light scalar mesons in the framework of the diquark chiral 
effective theory, which was introduced recently to describe exotic hadrons like pentaquarks. 
Using the parameters fixed previously by the mass of scalar mesons, we determine the axial 
couplings of diquarks with pions and kaons from the experimental data for the decay width 
of scalars in the dominant decay modes. The axial couplings turned out to be small, showing 
that the naive diquark model of Jaffe and Wilczek is reliable. Then, we calculate the decay 
width of scalar mesons for the sub dominant channel and find that our result is not far off 
from the experimental data. 

To conclude, in the diquark chiral effective theory, based on the spontaneously broken 
chiral symmetry and the diquark picture of exotic hadrons, we find the decay width of scalar 
mesons in good agreement of experimental data. It shows that the diquark chiral effective 
theory captures correct physics of exotic hadrons. 
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